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ABSTRACT 

The origin of dwarf spheroidal galaxies (dSphs) in the Local Group (LG) remains an enigma. The 
tidal stirring model posits that late- type, rotationally-supported dwarfs resembling present-day dwarf 
irregular (dirr) galaxies can transform into dSphs via interactions with Milky Way-sized hosts. Using 
collisionless A/'-body simulations, we investigate for the first time how tidal stirring depends on the 
dark matter (DM) density distribution in the central stellar region of the progenitor disky dwarf. 
Specifically, we explore various asymptotic inner slopes 7 of the dwarf DM density profiles {p oc 
as r ^ 0). For a given orbit inside the primary galaxy, rotationally-supported dwarfs embedded in 
DM halos with core-like density distributions (7 = 0.2) and mild density cusps (7 = 0.6) demonstrate 
a substantially enhanced likelihood and efficiency of transformation into dSphs compared to their 
counterparts with steeper DM density profiles (7 = 1). Such shallow DM distributions are akin to 
those of observed dirrs, highlighting tidal stirring as a plausible model for the origin of the morphology- 
density relation in the LG. When 7 < 1, a single pericentric passage can induce dSph formation and 
disky dwarfs on low-eccentricity or large-pericenter orbits are able to transform into dSphs; these new 
results allow the tidal stirring model to explain the existence of virtually all known dSphs across a 
wide range of distances from their hosts. A subset of rotationally-supported dwarfs initially embedded 
in DM halos with shallow density profiles are eventually disrupted by the primary galaxy; those that 
survive as dSphs are generally on orbits that are biased towards lower eccentricities and/or larger 
pericenters relative to those of typical cold dark matter (CDM) satellites. The latter could explain 
the rather peculiar orbits of several classic LG dSphs such as Fornax, Leo I, Tucana, and Cetus. We 
conclude that tidal stirring constitutes a prevalent evolutionary mechanism for shaping the nature of 
dwarf galaxies within the currently favored CDM cosmological paradigm. 

Subject headings: galaxies: dwarf - galaxies: Local Group - galaxies: kinematics and dynamics - 
galaxies: formation - methods: numerical 



1. INTRODUCTION 

The dwarf spheroidal (dSph) satellites in the Lo- 
cal Group (LG) are the faintest and most dark-matter 
(DM) dominated galaxies known (e.g., Mateo 1998; 
iTolstoy et al.l [20091 ) . These intriguing objects exhibit 
no appreciable gas components and tend to be clustered 
around the massive spiral galaxies Milky Way (MW) and 
M31. Furthermore, dSphs are characterized by pressure- 
supported, spheroidal stellar components (e.g., Mateo 
[T998 ) and a wide dive rsity of star formation histories 
(e.g-. lOrban et al.ll2008f ). Despite recent advances in our 
understanding of dSphs, a c onclusive model for their ori- 
gin has yet to emerge fe.g.. lKravtsQvl [2010). 

One class of models attempts to explain the forma- 
tion of dSphs via environmental proc esses, including tidal 
and ram pr essure stripping (e.g . . lEinasto et al.l 11974 
Faber fc LinI 119831: iMaver et al.l [200TI: Kravtsov et alJ 

[Kiimentowski et al. 
Epkas et al. 2011), 
(Kazantzidis et al. 
resonant stripping 
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(jP^Onghia et al.l I2Q09D. In th is context, the "tidal 
stirring" model of *M aver et al.l (^01) postulates that 
interactions with MW-sized host galaxies can trans- 
form late- type, rotationally-supported dwarfs resembling 
present-day dwarf irregular (dIrr) galaxies into pressure- 
supported, spheroidal stellar systems with the kinematic 
and structural pro perties of the classic LG dSphs (see 
iLokas et al.ll2012al for an extension of this model to the 
newly-discovered ultra- faint dSphs). 

Previous work on tidal stirring; has almost exclusively 
adopted cuspy, 'Navarro et al.' (1996, hereafter NFW) 
density profiles to model the DM halos of the progenitor 
disky dwar fs. Yet, rotation-curve modeling; of o bserved 
dIrrs fe.g.. IWeldrake et al.ll2003l : iQh et al.ll2011f ). which 
constitute the most plausible predecessors of dSphs ac- 
cording to the tidal stirring model, favors DM halos with 
significantly shallower inner density slopes. Such shallow 
DM density distributions are also supported by hydrody- 
nami cal cosmolog;ical s imulations of dwarf galaxy forma- 
tion (jGovernato et al.l [2010) and recent analytic models 
for the origin of co re-like DM density profiles in dwarfs 
(|Pontzen Govern ato 2012). 

Given that tidal stirring involves a combination of 
tidally-induced dynamical instabilities (e.g., bars) and 
impulsive tidal heating of the stellar component, the 
cuspiness of the dwarf DM density distribution in the 
central region where the stars reside should play a cru- 
cial role in this complex transformation mechanism. In- 



2 



KAZANTZIDIS, LOKAS, & MAYER 



deed, disky dwarfs with shallow DM density profiles 
are characterized by lower central densities and, corre- 
spondingly, longer internal dynamical times compared 
to their counterparts with steeper DM density distri- 
butions. Systems with longer dynamical times respond 
more impulsively to external ti dal perturbations an d suf- 
fer stronger tidal heating (e.g. jGnedin fc Ostrikerl FlQQQ). 
Rotationally-supported dwarf galaxies with progressively 
shallower DM density profiles may thus experience in- 
creasingly probable and efficient transformations into 
dSphs. Nonetheless, there has been no quantitative work 
aimed at investigating this qualitative expectation. 

Here we explore this issue via a series of tidal stirring 
simulations of disky dwarfs embedded in DM halos with 
different inner density distributions. Our results demon- 
strate that the likelihood and efficiency of transformation 
into a dSph are both enhanced substantially when the 
progenitor rotationally-supported dwarfs possess shallow 
DM density profiles, whose slopes in the central stellar 
regions are in agreement with those inferred from recent 
theoretical and observational efforts in dwarf galaxy for- 
mation. These findings further establish tidal stirring as 
a prevalent evolutionary mechanism for shaping the na- 
ture of dwarf galaxies within the currently favored cold 
dark matter (CDM) cosmological paradigm. 

2. METHODS 

A descrip ti on of t he adopted methodology is given in 
iLokas etld] ()2Q12af ). For completeness, we provide a 
summar y of our approach he re. We employed the tech- 
nique of IWidrow et al.l (|2QQ8f ) to construct numerical re- 
alizations of dwarf galaxies consisting of exponential stel- 
lar disks embedded in DM halos. The DM ha lo density 
profiles followed the general form (lLokasll2QQ2| ). 

p(r) = — (1) 

where p^, r^, and 7 denote the characteristic inner den- 
sity, the scale radius, and the asymptotic inner slope of 
the profile, respectively, ps depends on the epoch of halo 
formation, the present-day values of the cosmological pa- 
rameters, and 7 (throughout the paper we assume the 
concordance ACDM cosmogony and z = 0). 

To examine how the DM density distribution in the 
central stellar region of the progenitor disky dwarf af- 
fects tidal stirring, we varied 7 in three otherwise identi- 
cally initialized dwarf galaxies. Specifically, we adopted 
7 = 1 (corresponding to the NFW profile) and two shal- 
lower inner slopes, namely a mild density cusp (7 = 
0.6) and a nearly constant density core (7 = 0.2). 
These shallow power-law indices are well-motivated as 
they are akin to those of both ob served dirrs (e.g., 
IWeldrake et~all 120031 : lOh et al.l[20Tll ) and realistic dlrr- 
like syste ms formed in hydro dyn amical cosmological sim- 
ulations (jGovernato et al.ll2010 ^. 

Each dwarf galaxy comprised a DM halo with a virial 
mass of Mvir = lO^M© (corresponding to a virial radius 
of Tvir ~ 25.9 kpc) and a concentration parameter of 
c = Tvir/^s = 20 {rs ~ 1 .29 kpc), and was expone ntially 
truncated beyond rvir (|Kazantzidis et ahl l2004a[ ). All 
DM halos hosted an identical stellar disk whose mass, ra- 
dial scale- length, sech^ vertical scale- height, and central 
radial velocity dispersion were equal to Md = 0.02Mvir, 
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01 


125 


5 


25 


2.1 


02 


85 
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17 


1.3 


03 


250 
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50 


5.4 


04 


125 


10 


12.5 


1.8 


05 


125 


2.5 


50 


2.5 



Rd ~ 0.41 kpc , Zd = 0.2Rd, and ctrq = 10 kms~^, re- 
spectively (see I lKazantzidis et al.l l2011al for the motiva- 
tion behind these parameter values). We note that the 
value of Rd is derived assuming a dim ensionless halo spin 
parameter of A = 0.04 ( Mo et al.lll998f ). Being embedded 
in different DM density distributions, the resulting disks 
differed in their velocity dispersion profile and Toomre 
Q stability parameter (at R = 2.5Rd^ Q ^ 3.8, 3.3, and 
2.9 for 7 = 1, 0.6 and 0.2, respectively; direct numeri- 
cal simulations of the evolution of all dwarf galaxies in 
isolation for a period of 10 Gyr confirmed their stability 
against bar formation). 

Our choices of 7 above lead to DM density profiles, 
circular velocity profiles, and stellar binding energy dis- 
tributions that display markedly dissimilar shapes within 
5Rd ^ 2 kpc, namely the radius containing t he vast ma- 
jority 95%) of disk stars (see Figure 1 of lLokas et al.l 
l2012af ). Hence, although our modeling approach is cer- 
tainly not unique, it does satisfy the main require- 
ment of the present study by ensuring that the em- 
ployed rotationally-supported dwarfs exhibit substan- 
tially different initial properties in the central stellar re- 
gion. In particular, decreasing cusp slopes correspond 
to less concentrated mass and energy distributions, less 
steeply rising circular velocity profiles, and smaller max- 
imum circular velocities Fmax (^max ^ 19.6, 17.8, and 
16.7 kms~^ for 7 = 1, 0.6 and 0.2, respectively). There- 
fore, rotationally-supported dwarfs embedded in DM ha- 
los with shallow density profiles are expected to manifest 
rather different responses to tidal shocks and overall tidal 
evolutions compared to their counterparts with steeper 
DM density distributions. We shall quantify these differ- 
ences in the following section. 

The dwarf galaxy models contained Nh = 10^ DM 
and Nd = 5 x 10^ disk particles with a gravitational 
softening of Ch = 60 pc and Cd = 20 pc, respec- 
tively. Resolution was adequate to resolve all scales of 
interest. We assumed a single host represented by a 
self-gravitating , high-res olution numerical model of the 
MW (Kazant zidis et al.ll2Ql la). Each disky dwarf was 
placed on five bound orbits of varying sizes and ec- 
centricities inside the primary galaxy, for a total of 15 
experiment s. In detail, we emplo yed a subset of or- 
bits from Kazantzidis et al.j (|2011a[ ). whose parameters 
were motivated by both theoretical studies of the orbital 
distri butions of cosniologica l halos in MW-sized hosts 
(e.g., iDiemand et~all 120071 : iKlimentowski et al.l I2010D 
and observational w ork pertaining to LG dwarfs (e.g., 
iMcConnachig 120121 ). Tabled summarizes the adopted 
orbital parameters, including apocentric distances, rapo, 
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Fig. 1. — Evolution of mass in the DM (left panels) and stellar (right panels) components of the simulated disky dwarf galaxies as a 
function of time. Results are presented for orbits Ol (upper panels), 03 (middle panels), and 05 (lower panels). 7 denotes the asymptotic 
inner slope of the dwarf DM density distribution and the solid, dashed, and dotted lines correspond to 7 = 1.0, 7 = 0.6, and 7 = 0.2, 
respectively. DM (stellar) masses are computed within 0.7 kpc from the center of the dwarf (see text for details) and are normalized to the 
initial DM (stellar) mass enclosed within 0.7 kpc, Mq. In all panels, vertical lines specify pericentric passages and horizontal lines indicate 
mass loss of 30%, 60%, and 90% with respect to the initial value. For a given orbit inside the host galaxy, rotationally-supported dwarfs 
embedded in DM halos with shallow density profiles (7 < 1) suffer enhanced mass loss compared to their counterparts with steeper DM 
density distributions (7 = 1). 



eccentricities, rapo/^peri, pericentric distances, rperi, and 
orbital times, Torb- The dwarf galaxies, having consid- 
erably lower masses with respect to the host, do not 
experience dynamical friction (e.g., 'Colpi et al.l Il999f ) 
and, as a result, their orbital parameters remain vir- 
tually unchanged during the course of their evolution. 
In all cases, the dwarfs started at apocenter and were 
evolved for 10 Gyr inside the primary. The alignments 
between the internal angular momenta of the dwarfs, 
those of the primary disks and the orbital angular mo- 
menta w ere all always p rograde and equal to 45° (see 
[Kazantzidis et al.l l2011al regarding the implications of 
these choices). Lastly, the simulatio ns were pe rformed 
with the AT-body code PKDGRAV (;Stadelll200lf ). 

3. RESULTS 

We examined the global response of the disky dwarfs 
subject to the tidal field of the host galaxy via the time 
evolution of their masses, kinematics, and shapes. These 
properties were always computed within 0.7 kpc from the 
dwarf center, a radius which is approximately equal to 
the half-light radius, ri/2, of the initial disk. Adopting 
this well-defined, fixed scale facilitates meaningful com- 
parisons among different experiments and enables us to 
overcome the co mplications associated with determining 
tidal radii (e.g.. ' Read et al.|[2QQ6i) . 

The dwarf kinematics and shapes were quantified 



through the parameters Vrot/c'"* and c/a, where Vi-ot, 
(T*, c, and a denote the rotational velocity, the one- 
dimensional velocity dispersion, the minor axis, and the 
major axis of the stellar distribution, respectively. At 
each simulation output, the following procedure was re- 
peated. First, we determine the directions of the princi- 
pal axes and derive c/a using the moments of the inertia 
tensor. Subsequently, we introduce a spherical coordi- 
nate system (r, 6>, ^), where (j) is the angle coordinate 
in the xy plane oriented in such a way that the z-axis is 
along the minor axis of the stellar distribution. Lastly, 
we calculate Vrot around the minor axis Vrot = and 
compute the dispersions a^, cre^ and around the mean 
values; a*, which we adopt throughout the paper as a 
measure of the amount of random motions in the stars, 
is defined as a* = [{a^ ere ^ cr2)/3]i/^ 

Given the typical parameters associated with orbit 01, 
we use it as the basis for the comparison with the other 
experiments. In what follows, we only describe results 
pertaining to orbits 01, 03, and 05. This is for brevity 
and because the evolution of the rotationally-supported 
dwarfs on the tightest (02) and the most eccentric (04) 
of our orbits displays similar salient features to those of 
01. We emphasize that orbits 03 and 05 are charac- 
terized by larger pericenters and/or lower eccentricities 
compared to those of representative CDM satellites (e.g.. 
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TABLE 2 
Summary of Results 
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iDiemand et~ani2QQ7[) . 

Figure [1] shows the time evolution of both DM and 
stellar mass of the simulated disky dwarf galaxies as they 
orbit inside the primary. Apart from illustrating the con- 
tinuous stripping of the dwarf mass by the host tidal field, 
this figure highlights a number of additional generic fea- 
tures. First, masses decrease significantly at pericentric 
passages, where the intensity and variation of the time- 
dependent tidal forces are the strongest, and the tidal 
shocks occur. Between pericenters, masses remain re- 
markably constant, indicating that the dwarfs respond 
nearly adiabatically to the tidal field. Furthermore, the 
mass loss of the DM component greatly exceeds that of 
the stellar distribution, reflecting the more efficient tidal 
stripping of the extended DM halos of the dwarfs. 

Second, irrespective of orbit, disky dwarf galaxies em- 
bedded in DM halos with shallow density profiles (7 < 1) 
suffer augmented mass loss in both stellar and DM com- 
ponents compared to their counterp arts with steeper 
DM de nsity distri butions (7 = 1) (see iKazantzidis et alj 
l2Q04bl : [Peharrubia^^et^l] [2Q1QI for similar conclusions; 
these studies, however, performed simulations of the tidal 
evolution of dwarfs without stellar disks). In some cir- 
cumstances, the amount of mass loss is so dramatic that 
it leads to complete tidal disruption. Columns 4 and 5 
of Table [2] provide information regarding dwarf survival 
(as indicated by the presence of a self-bound entity at 
t = 10 Gyr) and list the time of disruption if it occurs, 
respectively. 

In Figure [21 we present the time evolution of the kine- 
matics and shapes of the simulated disky dwarfs as they 
orbit inside the host galaxy. Given that Vi-ot is larger for 
more concentrated mass distributions and that all dwarf 
models are initialized with the same value of ctro^ steeper 
cusp slopes result in greater initial values of Kot / c^* • We 
also note that the 7 = 0.2 dwarf exhibits Kot/c^* ^ 1 
initially. This notably low value of Vrot / cr^ for a disk is 
simply a consequence of our specific choice to compute 
the dwarf properties within 0.7 kpc from the center and 



does not interfere with the interpretation of the results. 

Injection of energy via tidal shocks at pericentric pas- 
sages increases the random motions of disk stars, causing 
their orbits to become more isotropic. Thus, Kot/c^* 
progressively decreases, signifying the transition from 
rotationally-supported to pressure-supported stellar sys- 
tems dominated by random motions. Simultaneously, 
the initially disky stellar distributions evolve into more 
spheroidal shapes as reflected in the continuous increase 
of c/a. Figure [2] shows that, regardless of orbit inside 
the primary, rotationally-supported dwarfs embedded in 
DM halos with shallow density profiles (7 < 1) demon- 
strate stronger evolution in their kinematics and shapes 
compared to their counterparts with steeper DM density 
distributions (7 = 1). 

It is important to stress that the strong tidal forces at 
pericenters typically trigger bar instabilities in the disks 
of the dwarfs. These tidally-induced bars play a vital 
role in the overall decrease of Vi-ot by transporting an- 
gular momentum outwards. As tidal stripping removes 
the outer parts of the dwarfs, the entire angular momen- 
tum content progressively decreases and the ability of 
the dwarf galaxies to be supported by rotation progres- 
sively diminishes. The seemingly unexpected increase 
of Vrot /cr^ observed in Figure [2] between some pericentric 
passages is due to the particular orientation of these bars 
at the moment of pericenter crossing and the strong tidal 
torques exerted on them by the primary galaxy, res ulting 
in an increase of Kot (see IKazantzidis et al.ll2011af ). Col- 
umn 6 of Table [2] indicates that bar formation is more 
probable in shallow DM density distributions (see also 
iMaver fc Wadslevl l2004D . 

Our primary goal is to assess the likehhood and effi- 
ciency of transformation into a dSph via tidal stirring. In 
accorda nce with observational studies of L G dwarf galax- 
ies fe.g. JMate3ll998l : [McConnachiell2012f ). we classify as 
bona fide dSphs only those simulated systems whose final 
states are characterized by Kot/c^* ^0-5 and c/a > 0.5. 
We underline that our theoretically derived values of 
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Fig. 2. — Evolution of stellar kinematics (left panels) and shape (right panels) of the simulated disky dwarf galaxies as a function of time. 
Results are presented for orbits Ol (upper panels), 03 (middle panels), and 05 (lower panels). Line types are as in Figure [T] Vrot /cr^ and 
c/a are computed within 0.7 kpc from the center of the dwarf (see text for details). Horizontal lines indicate the limiting values Vrot/cr* = 1, 
K-ot/c"* =0.5, and c/a = 0.5: simulated dwarfs whose final states are characterized by Vrot/cr* ^0.5 and c/a > 0.5 correspond to bona fide 
dSphs, while those with final values of 0.5 < Vrot/cr* < 1 and c/a > 0.5 are classified as "dSph?" (see text). For a given orbit inside the 
host, rotationally-supported dwarf galaxies embedded in DM halos with shallow density profiles (7 < 1) experience a stronger evolution in 
their shapes and kinematics and demonstrate a considerably enhanced likelihood and efficiency of transformation into dSph-like systems 
relative to their counterparts with steeper DM density distributions (7 = 1). 



Vrot /cr* and c/a are appropriate for m eaningful compar- 
isons w ith those of observed dSphs (seel Kazantzidis et al.l 
l2Qllal for a thorough discussion on this issue). In fact, 
the criterion Vrot/c^* ^ 0.5 for establishing dSph forma- 
tion in our experiments should be regarded as fairly con- 
servative. Indeed, we measure Kot around the minor 
axis of the stellar distribution, which corresponds to ob- 
serving the simulated dwarfs perfectly edge-on. Conse- 
quently, Vrot is nearly equivalent to the maximum rota- 
tional velocity and, thus, the values of Vrot/c^* that we 
quote throughout the paper are essentially the largest 
possible. Adopting a random line-of-sight would result 
in smaller Kot/c^* values, indicating even stronger and 
more complete transformations. 

Columns 7 and 8 of Table [2] list the final values of 
Vrot/o'^ and c/a, respectively. These values correspond 
to different timescales (from ~ 4 to 10 Gyr) that are 
determined by the presence of a well-defined, bound 
stellar component of ~ 0.7 kpc in size. Column 9 re- 
ports whether a dSph was produced according to the 
two criteria above. The notation "dSph?" in this col- 
umn is introduced to account for observed dSphs that 
exhibit 0.5 < V^rot/^* 5, 1 such as C etus ([Lewis et al.l 
I2007D. Tucana (iFraternah et al.l l 2009f ). and Andromeda 
II (|Ho et al.l l2012f ). Of the disky dwarfs initially em- 
bedded in cuspy DM halos (7 = 1), only those on 
high-eccentricity (rapo/'^peri ^ 5) and small-pericenter 



(^peri ^ 25 kpc) orbits are transfo rmed into objects 
with the properties of dS phs (see also'Maver et al.l l2001l : 
iKazantzidis et al.ll2011a[ ). Most importantly, the like- 
lihood of dSph formation is enhanced significantly for 
shallow DM density profiles (7 < 1); in such cases, 
rotationally-supported dwarfs on previously unfavorable 
low-eccentricity and/or large-pericenter orbits are able 
to transform into dSph-like systems. Figure [3] shows the 
surface density maps of the final stellar components of 
the dwarfs and visually confirms these conclusions. 

It is also worth emphasizing that, when 7 < 1, a dSph 
can be created without significant mass loss in the stel- 
lar component (see Figure [1]) because strong tidal shocks 
are not required to induce transformation for 7 = 0.6 
and 7 = 0.2. This finding has interesting implications 
as it suggests that the presence or absence of prominent 
tidal tails in present-day dSphs may not constitute a ro- 
bust observational constraint on the tidal stirring model 
without prior knowledg e of the inner DM density distri- 
butions of such objects (jMaver et al.ll2002l ). 

Inspection of Columns 6 and 9 of Table [2] illus- 
trates that dSph formation only occurs in conjunc- 
tion with the development of bar instabilities, high- 
lighting a strong association between bar formation 
and transformation into a d Sph via tidal stirring (see 
also, e.g., [ Maver et al.l 120011 : iKlimentowski et al.l l200l 
IKazantzidis et al. 201 lal ). Lastly, Column 10 of Table E] 
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Fig. 3. — Surface number density distribution of stars in the simulated dwarfs. Results are presented for orbits Ol (upper panels), 03 
(middle panels), and 05 (lower panels). All stellar components are depicted at t = 10 Gyr, except for that in the upper rightmost panel 
(Ol; 7 = 0.2) which is shown at t = 6.35 Gyr (corresponding to the timescale of the last apocenter before the disruption of the dwarf 
galaxy). The asymptotic inner slope of the dwarf DM density profile, 7, is indicated in each panel. The images show the most non-spherical 
views along the intermediate axis of the stellar distribution, y (where x and z denote the major and minor axes, respectively). The stars 
are binned into 0.2 kpc x 0.2 kpc fields perpendicular to the line-of-sight. The contours correspond to the number of stars within each 
bin and are equally spaced by 0.2 in log A^. The innermost contours are in the range log A" = 3.6 — 4.6 depending on the case. 



lists the time elapsed from the beginning of the simula- 
tion until dSph formation occurs (the number of corre- 
sponding pericentric passages is included in parentheses). 
The entries in this column indicate, that similar to the 
likelihood, the efficiency of transformation into a dSph 
is considerably increased for disky dwarfs embedded in 
DM halos with shallow density profiles. 

4. DISCUSSION 

The present study is the first to elucidate the effect 
of the asymptotic inner slope 7 of the dwarf DM den- 
sity profile (p oc as r ^ 0) on the tidal evolution 
of dwarf galaxies consisting of baryonic components in 
the form of stellar disks. We have shown that, regard- 
less of orbit inside the primary, rotationally-supported 
dwarfs embedded in DM halos with core-like density dis- 
tributions (7 = 0.2) and mild density cusps (7 = 0.6) 
demonstrate a substantially enhanced likelihood and ef- 
ficiency of transformation into dSphs relative to their 
counterparts with NFW-like steeper DM density profiles 



(7 = 1). Such shallow DM distributions are akin to those 
inferred from the mass modeling of observed LG dirrs 
(e.g., Weldrake et al. 2003) and an alogous galaxies out- 
side of the LG fe.g-. lOh et al.ll2011[ ). This fact highlights 
tidal stirring as a plausible mechanis m for the orig in of 
the morphology-density relation (e.g.. iMateol Il998f ). an 
essential constraint that any model for the LG must sat- 
isfy. 

The following order-of-magnitude calculations offer in- 
sight into our numerical results. In the impulse approxi- 
mation, the energy injected at each pericentric passage is 
given by /S.E oc M^^^^M R'^V~^ , where Mhost is the mass 
of the host enclosed within rped, M denotes the mass of 
the dwarf within a characteristic radius i?, and Vi-ei is 
the relative velocity of the two galaxies at the pericen- 
ter of the orbit (e.g., iBinnev fc Tremaiii3l2008f ). For a 
specific set of initial orbital parameters, Mhost and Vrei 
are fairly similar for different 7. Thus, AE oc MR^. 
Moreover, by virtue of the virial theorem, the energy 
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content of the dwarf scales as £^ cx jR. Hence, the 
fractional increase in energy caused by the tidal shocks is 
/S.E / E (X /M . At a given distance R from the center 
of our initial disky dwarf galaxies, decreasing cusp slopes 
correspond to smaller M, and thus to larger /S.E/ E. This 
explains why the 7 = 0.6 and 7 = 0.2 rotationally- 
supported dwarfs experience stronger tidal shocks and 
augmented mass loss relative to their 7 = 1 counterparts, 
leading to their enhanced morphological transformation 
into dSphs. Considering adiabatic corrections to the 
energy change predicted by the impulse approximation 
(e.g., iGnedin fc Ostriker 1999i) would only reinforce this 
conclusion. Indeed, owing to the fact that such correc- 
tions are inversely proportional to a power of the stellar 
orbital frequency uo and uj is higher at a given radius for 
more concentrated mass distributions, decreasing cusp 
slopes would correspond to even larger /S.E / E compared 
to those predicted by the impulse approximation. 

The results presented in this paper demonstrate that 
the orbital parameters and the DM density distribu- 
tions of the progenitor rotationally-supported dwarfs can 
independently determine the final properties of dSphs 
formed via tidal stirring. Therefore, the fact that Fornax 
and Draco, for example, have roughly similar masses at 
prese nt, as inferred frorn their stellar velocity dispersions 
(e.g.. iKazantzidis et al.ll2004bi ). but differ by about two 
orders of magnitude in luminosity, could be explained 
in two ways. One possibility is that the predecessors of 
these two dwarf galaxies acquired very different DM den- 
sity profiles, for reasons related to their formation history 
and not to the environment. Alternatively, Fornax and 
Draco may have originated from systems that had com- 
parable DM and baryonic mass distributions, but dis- 
played dissimilar tidal evolutions and evolved differently 
because they entered their host galaxy on different orbits, 
perhaps due to a different infall epoch onto the primary. 
Other effects not included here, such as stripping due 
to ram pressure and external radiation fields at high red- 
shift, may also affect the disky dwarfs differently depend- 
ing on their orbits (Mayer et al. 2006, 2007) and internal 
mass distributions, increasing further the scatter in the 
properties of the resulting dSphs. 

The increased mass loss and rate of disruption expe- 
rienced by dwarfs embedded in DM halos with shallow 
density profiles have important implica tions for alleviat- 
ing the ^'miss in^ satellites problem" (Moo re et al.l 119991 : 
iKlypin et al. 1^1999). These findings also suggest that the 
dwarf luminosity and mass functions must be shifted 
lower relative to those of ACDM simulations that do 
not take into account baryonic effects (such as outfiows 
triggered by supernovae explosions; see iGovernato et al.l 
[2010) which can give rise to fiattening of DM cusps. Con- 
sequently, the recent claims regarding the agreement be- 
tween theory and observations of the luminosity function 
and the radial d istribution of satelli tes within MW-sized 
host halos (e.g.. iMaccio "eTallllnoh may need revision. 

Our investigation establishes tidal stirring as an even 
more prevalent transformation mechanism than previ- 
ously considered. Indeed, earlier work adopting cuspy 
DM density profiles for the progenitor disky dwarfs has 
indicated that only orbits that are characterized by high 
eccentricities (rapo/^peri ^ 5), typical of CDM struc- 
ture formation models (e.g.. [pTemand et al.ll2007l ). and 
small pericentric distances (rperi ^ 25 kpc) are capable 



of tra nsforming rotationally-supported dwarfs int o dSphs 
fe.g.. lMaver et al.ll200ll : IKazantzidis et al.ll2011al ). For a 
transformation to occur, the same studies have also con- 
cluded that Torb should be short enough to allow the 
disky dwarf galaxies to complete at least two pericen- 
tric passages inside their hosts. These two requirements 
are probably not met for a number of LG dSphs, a fact 
that questioned the widespread applicability of the tidal 
stirring model. 

For example, the distant dSphs Leo I, Cetus, and Tu- 
cana (residing in the outskirts of the LG at several hun- 
dred kpc from the MW and M31; e.g., (C aputo e t al.l 
119991 : iMcConnachie et al.l[2QQ5l : iSaviane et al. .1996) are 
likely moving on very wide orbits associated with ex- 
tremely long Torb- Moreover, proper-motion measure- 
ments have indi cated a fair ly large pericenter for Leo I 
(^peri ^ 60 kpc; 'So hn et an '2012) and a low-eccentricity 
orbit with a similarly great p ericenter for the dSph For- 
nax (e.g., iPiatek et al.ll2007l ). Our results suggest that 
the origin of all these dSphs can be accommodated within 
the framework of the tidal stirring model (for alterna- 
tive formation scenarios of the distant LG dSphs and 
Fornax, see, e.g., Kravtsov et al. 2004; Sales et al. 20071: 
i Ludlow et al. 2009: IKazantzidis et al.l I2011bl and e.g., 
lYozin fc BekkLi2012 l. respectively). Indeed, rotationally- 
supported dwarf galaxies on low-eccentricity and/or 
large-pericenter orbits may be efficiently transformed 
into dSph-type objects, provided that their DM halos 
possess shallow density profiles (see Table [2j). In fact, 
assuming core-like DM distributions (7 = 0.2) and small 
^peri, even a single pericentric passage can induce trans- 
formation into a dSph. We note that according to recent 
studies, Leo I has pas sed through the pericenter of its or- 
bit once (Sohn et al. 2012) and the same conclusion may 
also apply to Cetus and Tucana (Teyssier et al. 2012). 

The present work confirms the previously reported 
strong association between the development of bar in- 
stabilities in the disks of the progeni tor late- type dwarfs 
and the formation of dSp hs (e.g., *M aver eFal 120011 : 
iKlimentowski et al.l 120091 : IKazantzidis et al.l \2mT^ . 
Bar-like structures should thus be common among the 
less evolved dwarf galaxies in the LG as the bar stage 
constitutes one of the longest phases in the transforma- 
tion process. Recently, the first steps toward validating 
this fundamental prediction of tidal stirring have been 
taken with the detection of bars in a number of MW 
satellites, including Ursa Minor, Sagittarius, and Carina 
(|Lokas et al.ll2012b[ ). Nonetheless, with the notable ex- 
ception of the Large Magellanic Cloud, the number of 
irrefutable bar-like distortions observed in LG dwarfs 



is sti ll relatively low (see, however, IKlim entowski et al.l 
I2OO9I for a discussion pertaining to the intrinsic difficul- 
ties in identifying bars in the dwarf galaxies of the LG). 

Our simulations also indicate that none of the 7 = 1 
disky dwarfs are disrupted by the tidal field of the pri- 
mary galaxy; of these, only systems on high-eccentricity 



(^apo/'^] 



apo/ ' peri ^ 



5) and small-pericenter (rperi ^ 25 kpc) 
orbits are transformed into objects with the properties 
of dSphs. Interestingly, when 7 < 1, the dSph-like sys- 
tems that survive are generally on orbits with lower ec- 
centricities and/or larger pericen ters compared to thos e 
of typical CDM satellites (e.g., iDiemand et al.l I2007D . 
This novel finding provides a natural explanation for the 
rather peculiar orbits of several classic LG dSphs such as 
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Fornax, Leo I, Tucana, and Cetus. 

The degree of transformation into a dSph depends 
on both the number and the strength of tidal shocks, 
which are determin ed by Tprb and rneri , respe ctively (e.g., 
[Mayer et al. 2001J; iKazantzidis et al.l [2011a). Specifi- 
cally, short orbital times and small pericentric distances, 
corresponding to orbits associated with a large number of 
strong tidal shocks, produce the strongest and most com- 
plete transformations. Assuming that the distant dSphs 
in the LG are products of tidal stirring and that their 
progenitors have no bias towards significantly shallower 
than average inner DM density distributions, we predict 
that these dwarfs should exhibit higher values of Kot / 
relative to those of dSphs located closer to their hosts. 
Although tentative evidence of intrins ic rotation exists in 
Leo I (e.g., Koch et al. 20 03), Cetus ([Lewis et al. 2007), 
and Tucana (Fraternali et al. 2009), future conclusive 
measurements of kinematics in these dSphs will serve to 
validate (or falsify) this prediction. 

Lastly, our approach neglects the effects of gas dynam- 
ics, star formation, and chemical evolution. Definitive 



conclusions regarding the efficiency of the transforma- 
tion process clearly require more sophisticated theoreti- 
cal modeling, capable of capturing the turbulent, multi- 
phase structure of the interstellar medium in dwarf galax- 
ies. Investigating the tidal stirring of realistic dirr galax- 
ies, such as those formed self-consisten tly in the cosmo- 
logical hydrodynamical simulations of [Governato et al.[ 
(|20T0[ ). with recipes of radiative cooling, star formation, 
and supernovae feedback would thus be particularly valu- 
able in this direction. 
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